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Abstract. Bayesian Blocks is a new time series algorithm for detecting localized 
structures (spikes or shots), revealing pulse shapes, and generally characterizing in- 
^-v \ tensity variations. It maps raw counting data into a maximum likelihood piecewise 

£»\j ■ constant representation of the underlying signal. This bin-free method imposes no 

lower limit on measurable time scales. Applied to BATSE TTE data, it reveals the 
shortest know burst structure - a spike superimposed on the main burst in GRB 920229 
= Trigger 1453, with rise and decay timescales ~ fewlOO/isec. 
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While there is phenomenological understanding of some general features of the in- 
credibly varied GRB time profiles [10,16,6,1,3,17], the underlying physical processes 
are unknown. We are analyzing BATSE data with a variety of techniques aimed 
at extracting physically relevant information without regard to physical models. 

The BATSE time-tagged event, or TTE, data have the highest time resolution 
(but are perforce limited to short bursts). We have developed a bin- free time 
domain method, based on a rigorous statistical treatment of the Poisson photon 
detection process, to characterize the underlying signal in spite of the large sta- 
tistical fluctuations inherent in photon counting. Called Bayesian Blocks (BB) 
on account of the nature of its statistical foundation, the technique starts from 
raw photon counts and produces the most likely piecewise constant representation 
of brightness as a function of time. In this paper we treat only TTE data, but 
elsewhere the same procedure is developed for binned and TTS data [12]. 



THE ANALYSIS METHOD: BAYESIAN BLOCKS 

The idea rests on the simplicity of mathematical models of two processes: a 
source of constant brightness over a given time interval (this model is here denoted 
Mi), and the detection of an individual photon. These simplicities lead to an 
exact formula for the model likelihood for a given array of photon arrival times. 
Marginalization of nuisance parameters with physically derived priors [12,13] yields 
the global likelihood 

where N is the number of photons detected in the interval, M is the length of the 
interval in the units in which arrival times are quantized, and r(iV + 1) = N\ is the 
gamma function. When compared to the analogous global likelihood of a second 
model, this quantity measures the probability of the constant source model. 

Indeed, the other model we are interested in, M2, is that in which the observed 
interval is divided into two subintervals, separated by a change point, each of which 
has a constant Poisson rate. The global likelihood of this compound model is 

L(M 2 \D TTE ) = £ 0(JVi, M 1 )(t>(N 2 , M 2 ) (2) 

where the pairs N, M are the number of photons in, and length of, each subinterval, 
and the sum over n cp marginalizes the change point location. 

Taking p to be the ratio of the prior probabilities of segmented and unsegmented 
models, the odds ratio 

L(M 2 \D TTE ) 

(J = P-, : (3) 

r £ (Mi I D TTE ) 

gives the relative posterior probabilities of the two models (for the same data). The 
choice p = 1 expresses lack of prior knowledge of the structure. If the odds ratio 
in (3) favors segmentation, the location of the optimal change point is found by 
maximizing the summand in the eq. (2). This is easy, because the summand can 
be written as a function of only one variable - the location of the photon at the 
change point. In turn the maximum likelihood Poisson rates are found to be j^ for 
each subinterval. The full procedure to form a block representation of the signal 
consists of iterative application of the decision making step just described, until 
each subinterval is judged better off left undivided. 

We are using this algorithm [11] to detect and characterize pulse structures known 
to make up the time-profiles of many 7-ray bursts [10]. The BB decomposition 
has proved to be an excellent, completely automatic, way to determine approxi- 
mate widths, locations, and amplitudes of pulses, without invoking parametric or 
other explicit pulse-shape models. These are fed to a nonlinear optimizer which 
deconvolves, and estimates parametric models of, overlapping pulses. Moreover the 
information contained in the blocks is itself useful, especially in cases where the 
overlap is not too severe. 



A MICROBURST IN GRB 920229 = TRIGGER 1453 



Trigger 1453 nicely displays the time resolution and dynamic range of the BB 
procedure. From the time of the First BATSE Catalog [4] it has been commonly 
known that this burst contains a spike of duration <m 10 ms near the end of a 
relatively flat burst of duration ps 0.192 sec. 
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FIGURE 1. Bayesian Block decomposition for Trigger 1453, all energy channels. Count rate 
(in units of 1000 c/sec) vs time (sec). BB: thick lines; binned TTE data (for display only): thin 
lines. The microburst peak at 0.191s sec is offscale. The last BB shown extends to the end of the 
data at 1.42 sec, and is a good estimate of the background, (a) p = 1; (b) p = .01. 



With no special treatment this spike showed up as the shortest in the distribution 
of block widths. Figure 1 shows the BB decomposition of the TTE data. The the 
top panel, for p = 1, shows tall thin blocks at several places on the broad plateau. 
These undoubtedly spurious whiskers are probably signs of occasional errors in the 
segmentation decision, based on O > 1 vs O <= 1 in eq. 3; for O ~ 1 this coin 
flip should go wrong about half the time. Taking p < 1 wards off these errors, as 



seen in the bottom panel; real structures, including the (offscale) spike, typically 
correspond to astronomically large odds ratios and are therefore barely affected. 

Figure 2 includes individual energy channels and expands the region of the spike 
- confirming its reality, except in the low S/N channel 4. This is the only short 
spike in our BB analysis of 317 short bursts, in agreement with a search for spikes 
in the range 10/is — 10ms in TTE data for 20 bright bursts [15]. 
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FIGURE 2. The microburst as in Fig 1; p = 1. (a) sum of energy channels; (b-d) channels 1-4. 



OTHER WORK 



A more detailed study of this burst will address - using a sensitive Bayesian 
method [18] for assessing the likelihood that data histograms are from the same 
distribution - whether the detector (or energy) distributions of spike photons are 
different from those of the main burst. This seems unlikely, but would argue that 
the spike, in a different part of the sky, is different from ordinary bursts. 



We are investigating replacement of the ad hoc iteration used here with the more 
rigorous approach of multiple change point detection [2]. While this may lead 
to improvements, numerical experiments with real and synthetic data show our 
implementation of BB to be accurate and efficient. Other relevant work includes 
a negative search for microsecond flares [14] , a study of fundamental instrumental 
limits on detection of short spikes [9] , a Bayesian method to detect bursts in binned 
count data [7], and wavelet methods [5]. 

Thanks to Tom Loredo, for many suggestions, and NASA's Astrophysics Data 
Program. 
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